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Recent studies of azo-dye doped liquid crystal elastomers show a strong photomechanical response. 
We report on models that predict experimental results that suggest photothermal heating is the 
dominant mechanism in a planar constrained geometry. We compare our models with experiments to 
determine key material parameters, which are used to predict the dynamical response as a function 
of intensity. We show that a local strain from photothermal heating and a nonlocal strain from 
thermal diffusion is responsible for the observed length changes over time. This work both elucidates 
the fundamental mechanisms and provides input for the design of photomechanical optical devices, 
which have been shown to have the appropriate properties for making smart materials. 
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I. INTRODUCTION 

There has been a growing interest in photomechani- 
cal behavior of organic materials over the last decade. [TL 
Q Recently, an experimental study has determined the 
mechanisms of photo- induced deformations of liquid crys- 
tal elastomers in a photomechanical optical device (POD) 
geometry. Q This work was motivated in part by the cre- 
ation of parallel beam PODs made of high-concentration 
azo-dye doped LCEs that were cascaded together in series 
to show a proof-of-concept network of sensor/actuator 
photomechanical devices. [7[ 

The observation of the photomechanical effect dates 
back to Alexander Graham Bell, who showed that a voice 
could be transmitted over a beam of light, [sl] A century 
later, Uchino et al demonstrated the "Uchino walker." 0- 
IT]| This device implimented ceramics that constrict when 
exposed to light, and relax in the dark (reversible pho- 
tostriction). The mechanism of ceramic photostriction 
can be explained by light-induced charge diffusion fol- 
lowed by a piezo-electric effect, which is induced by the 
resulting internal electric field. 

With breakthroughs in the technology for making 
polymer optical fibers, the first all-optical photome- 
chanical device was constructed that acted as a posi- 
tion stabilizer. This device encompasses five device 
classes: information processing, logic, transmission, sens- 
ing, and actuation, [ij] The next year, a small-scale, 
2.5 cm, photomechanical optical device was made using 
an azobenzene-dye-doped multi-mode polyme r op tical 
fiber with frayed ends that act as retroreffectors.fl^ This 
device was later used to demonstrate all-optical modula- 



tion of a probe beam by introducing a pump beam that 
changed the fiber length, [l^ A decade later, a fiber can- 
tilever was observed by the use of off-center propagation 
of light throu gh a disperse red 1 dye-doped polymer op- 
tical fiber. [11,13 

Recently, nematic liquid crystal elastomers (LCE) have 
been reported to undergo large deformations via photo- 
isomerization. (TtI [18| A high-concentration azobenzene- 
dye-doped nematic LCE was shown to swim over the 
surface of water when illuminated by a transverse light 
source, [l^ Other studies of nematic LCEs have investi- 
gated the nonlinear behavior of a LCE cantilever. (20| 

Motivated by the recent experimental observations of 
a nematic LCE in a POD geometry, this work seeks to 
build an understanding of the mechanisms. We present 
numerical modeling results of the mechanisms and com- 
pare the results with experimental observations of the 
time-dependent changes in length of a LCE in response 
to light. Understanding these mechanisms is the first 
step in the creation of ultra-smart materials, plj 



II. PHOTOMECHANICAL OPTICAL DEVICE 
CONFIGURATION 

We develop a theory that describes the mechanisms of 
photo-induced deformations of liquid crystal elastomers 
in the surface-constrained geometry as described below. 
During these experiments, a LCE is pressed between two 
parallel glass substrates with a small force, which is ap- 
plied to achieve a visible "wet spot" at the LCE-glass 
interface. Silver is deposited on the inner surfaces of the 
glass - away from the position of the LCE. The pump 
beam is steered into the glass substrate and incident upon 
one end of the LCE. The probe beam is directed through 
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FIG. 1. (Color online) Diagram of the photomechanical opti- 
cal device (POD) used to measure the LCE's length change 
by detecting changes in the interference pattern of the probe 
beam. 



the interferometer region of the POD, i.e. where the sil- 
ver is deposited. 

The uniaxial nematic LCE is oriented between the two 
glass substrates so that the director orientation is normal 
to the LCE-glass interface. The LCE acts as one of the 
three stabilizing points that define the planar surface of 
one substrate. The other two points are defined by the 
tips of adjustable thumb screws. These were used to 
bring the two substrates into parallel alignment. Figure [1] 
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FIG. 2. The chemical structures of the silicon backbone, tri- 
functional crosslinker, mesogenic sidechain, and disperse or- 
ange 3 dopant chromophore that are used to construct the 
dye-doped LCEs. 




LCE 



Xq X2 X3 X4 



FIG. 3. (Color online) A liquid crystal elastomer illuminated 
by a 488nm wavelength laser, where the direction of light 
propagation is parallel to the director orientation. The gra- 
dient represents laser light absorption. 

shows a perspective diagram of the POD used in previous 
experiments. [6| 

The LCE material within the POD consists of a silicon 
backbone, crosslinker, and mesogenic sidechain whose 
structures are shown in Figure [2] The sidechain is con- 
nected to the backbone by the crosslinker, and there 
are roughly 10 silicon backbone segments per crosslinker. 
Disperse orange 3 (DOS) is a photo-isomerable dye that 
is dissolved in the LCE to act as an absorber of the 
pump beam. These absorbers can both change shape 
and transfer heat to the surrounding mesogens to reduce 
the macroscopic orientational order. [2^ 

Our models treat the high concentration limit, i.e. 
approximately 0.1% by weight DOS-doped LCEs. The 
LCEs are pumped with light at a frequency of 488 nm, 
near the resonant absorption peak.Q At this concentra- 
tion and wavelength, fi{Q = 0.7) w 0.2 /im~^, where /i is 
the Beer-Lambert coefficient (defined in the next section) 
and Q is the scalar uniaxial order parameter. 

III. THEORY 

It is well known that a continuous differential equation 
can be approximated by finite differences provided that 
the finite pieces are smaller than the smallest character- 
istic dimension. Figure |3] shows a schematic diagram of 




FIG. 4. (Color online) Sections of an LCE (top) before and 
(bottom) after turning on the laser. The sections decrease by 
a length, ALi, at position Xi. 
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an illuminated LCE in one dimension. Each mathemat- 
ically discrete section of the LCE experiences a length 
change that is much smaller than the length of each sec- 
tion, which allows for an engineering strain approxima- 
tion. A diagram is shown in Figure |4l 

A. Photo-isomerization induced strain 

The LCE's length at coordinate Xi changes by an 
amount ALj depending on the local light intensity and 
temperature, leading to a strain ALi/Ax. For the un- 
elongated elastomer, the Beer-Lambert law gives 

/ - /oe— (1) 

where / is the transmitted intensity, Iq the incident in- 
tensity, e the absorption coefficient, A the illuminated 
area, c the concentration, and I the path length. Fig- 
ure [3] shows that the intensity is maximum at the surface 

X = Xq. 

We first consider the effects of photo-isomerization, 
which is a local mechanism that results in a length change 
of element i in proportion to the intensity in element 
i. It is well known that the conformation of an azo-dye 
can change from trans to cis by photo-isomerization. [23j 
When azobenzene molecules are dissolved into a LCE, 
the process of photo-isomerization can cause a reduction 
in the long-range orientational order. [24] 

The probability that an azo-dye in the LCE under- 
gos trans-cis photo-isomerization is proportional to the 
population of trans isomers oriented aloiig_the light's po- 
larization axis and the light's intensity. [l8| We define ^ 
as the probability per unit of intensity per unit of time 
that a trans isomer will absorb light and be converted 
to a cis isomer. The temperature-dependent decay rate 
of the cis population, /3, is assumed to be independent 
of intensity. As long as the temperature increase in the 
elastomer is small, which it is in our experiment, this will 
be a good approximation. 

The photothermal heating and photo-isomerization 
mechanisms are coupled because the heating rate de- 
pends on the trans population, which depends on the 
rate of photo-isomerization. If the LCE is highly-doped, 
so that the 1/e absorption length is much smaller than 
the length of the sample, only the region near the sur- 
face of the sample will be illuminated. Thus, photo- 
isomerization-induced strain will be limited to a negli- 
gibly small part of the sample, where even a large local 
strain would produce only a small change in the total 
length of the LCE. 

Since the absorption coefficient decreases as trans iso- 
mers are converted to cis isomers, the laser penetrates 
deeper into the sample over time. The equilibrium pene- 
tration depth will be determined by the equilibrium trans 
population, which occurs when the photo-isomerization 
rate balances the recovery rate, and the difference in op- 
tical absorption coefficient of the cis and trans isomer. 
We find that this penetration depth, the 1/e absorption 



length, is approximately 5 ^m for the range of intensities 
used in our experiments at a wavelength of 488 nm. As 
such, we can treat the illuminated region as a thin heat 
source that contributes little to the total length change. 

The temperature of the rest of the LCE will increase 
due to thermal diffusion of heat from the illuminated re- 
gion, which leads to a reduction in the long-range ori- 
entational order parameter away from the optical heat 
source, and thus a change in the length. By virtue of the 
fact that the LCE is much longer than the illuminated 
region, even a small strain in the dark region of a highly 
azobenzene-doped LCE can lead to a length change that 
is much larger than in the illuminated region. 

Thus, the physical picture of the process is as follows. 
The light is absorbed in a thin layer of the LCE near 
its surface. The heat generated in this thin region fiows 
through the sample, causing a decrease in the orienta- 
tional order of the liquid crystalline molecules due to the 
resulting change in temperature. The order reduction 
can be explained using statistical models of nematic liq- 
uid crystals such as the method described by Maier and 
Saupe.fli-Il^ 

The change in orientational order causes a length 
change. Note that the population of cis isomers of 
the azobenzene dopants slightly increases with increased 
temperature, which at extremely high concentrations of 
photonematogen dopants, could also contribute to the 
decrease of orientational order of the liquid crystal, and 
thus also contribute to length contraction. We lump 
these two processes together and call it simply the ther- 
mal mechanism. The sample reaches thermal equilibrium 
when the heat fiowing out of the LCE into the glass sub- 
strates and lost due to convection at the glass-air inter- 
face is in balance with the heat generated in the thin 
surface layer. 

Under the assumption that the temperature change 
is small enough in the illuminated region to invoke the 
decoupled mechanism approximation, the population dy- 
namics are modeled by light-induced depletion of trans 
isomers and constant-rate recovery. The rate of change in 
the population fraction of trans isomers, N, throughout 
the LCE via photo-isomerization is 

dN (r f) 

]^ ' ' = -^N {x, t) I (x, t) + /3 (iVeq - N {x, t)) , (2) 

where A^eq is the equilibrium population fraction of trans 
isomers under ambient conditions and / (x, t) is the light 
intensity as a function of depth and time. While an in- 
crease in the temperature can both cause conversion from 
cis to trans as well as trans to cis isomers, in equilibrium, 
the population of the higher-energy species will increase 
with temperature. 

The intensity as a function of depth at the entry surface 
of the LCE is time-independent, that is, l{0,t) = Iq. 
Substituting the Beer-Lambert law into Equation [5] gives 
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the nonlinear rate equation 

d [ r ' 

—N(x,t)^-^N(x,t)Ioexp -/i / N(x',t)dx' 
at V Jo . 

+PiN,q~N{x,t)) (3) 

where fi — ecA is the absorption coefficient multiplied by 
the linear molecular density. 

Equation [3] is a generalization of the two state photo- 
isomerization rate equation in the small absorption limit 
given by Bian et al for photo- reorientation, 'l6| 



and 



dN 

^ = -e/A^ + /3(l-2iV), 



(4) 



which applies to the case when the material is much 
shorter than the absorption length. The 2N term ac- 
counts for molecules being equally distributed along two 
orthogonal directions in a plane perpendicular to the 
beam polarization axis. 

The strain induced by photo-isomerization, tXp, is 



(5) 



where ALp is the change in length caused by photo- 
isomerization and Lq is the initial length. 

The order parameter for a LCE is a function of the iso- 
mer concentration and temperature. The equation gov- 
erning the order parameter as a function of trans isomers 
is calculated self consistently, [l^l With a small cis iso- 
mer population at a given depth within the LCE, the 
uniaxial scalar order parameter, Q, can be approximated 
as a linear function of the dopant concentration, N , or 
AQ oc AA^, where A denotes the change in a variable. 
Also, the curve defining the strain as a function of change 
in the order parameter will be locally linear over small 
changes in order. Thus, the strain can then be expressed 
as at (x, t) oc AQ. This approximation allows the photo- 
isomerization strain to be expressed as 

ap{x,t) = ~b{N,^~N), (6) 

where h is the proportionality constant of photo- 
isomerization-induced strain. 



B. Temperature induced strain 



T(a;,0) = <?(x), 



(10) 



where -fCLCE is the LCE's thermal diffusivity, fcLCE is 
the LCE's thermal conductivity, T is the temperature, 
Hs {x,t) is the rate of temperature increase due to ab- 
sorption of light in the material, / describes the temper- 
ature gradient as a function of time at each boundary, rix 
is the unit vector component along the x direction, and 
g is the initial temperature profile of the material for a 
LCE of length Lq. 

Solutions to the boundary value problem depend on 
the experimental constraints and material parameters. 
For a steady initial state with no heat source, the tem- 
perature profile is constant. Thus, g = Tq (room tem- 
perature). Since the change in temperature is of inter- 
est rather than the absolute temperature, we subtract 
the ambient temperature from all temperatures, thereby 
defining Tq = 0. The light is turned on at t = to a con- 
stant intensity 7o> and remains so indefinitely. This step 
function intensity will result in transient behavior until 
the system is in equilibrium. A fit of the transient be- 
havior to the model that follows will allow the material's 
photomechanical parameters to be determined. 

The absorption of light from the laser acts as the heat 
source, so the amount of heating at a point in the LCE 
is proportional to the light intensity at that point. The 
boundary function / will be described by contact conduc- 
tion at the LCE-glass interface and the contact conduc- 
tance coefficient will be approximated as constant over 
small temperature increases. Therefore, the boundary 
conditions may be written as 



d 
drix 
Ci 

fcLi 



CE 



Tlce {x, t) 

X- 

[rLCE(0,t) 



Tglass (0,i)] 



(11) 



The other photomechanical mechanism that we con- 
sider is photothermal heating. The temperature dis- 
tribution throughout the LCE is described by the heat 
equation. The one-dimensional heat equation with a heat 
source and conduction-type boundary conditions is given 
by 



lT{x,t)-K^CE^T{x,t) 

khCE -^T{x,t) 
arLT 



fcLCE 3 — T{x,t) 



Hs{x,t), (7) 
= /i it) , (8) 
= /2 (t) , (9) 



and 



3 — Tlce (x, t) 
dn^ 



x=Lo 



^-LCE 



[TlcE (io, t) - Tglass (io, t)] 



(12) 



Assuming that both boundaries have the same contact 
conductance coefficient at the LCE-glass interface (Ci = 
C2 = C), the heat equation with the above boundary 
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conditions and the initial condition become 

-J^ThCE [X,t) - i^LCE^^TLCE [X,t) , 



3 — Tlce (a;, t) 



x=0 



c 



clce 



[TlCE (0,t)-Tgiass(0,t)], 



3 — Tlce (a;, t) 



x=Lo 



c 



^LCE 

Tlce (a:,0) = 



(13) 



(14) 



(15) 
(16) 



The light energy absorbed at any point in the mate- 
rial is proportional to the intensity at that specific point 
at that time and is also proportional to the number of 
absorbers oriented along the polarization of the light at 
that point in the LCE. We define ^ to be the constant of 
proportionality that relates the heating rate, Hs, to the 
intensity and trans population fraction, 



Hs{x,t) = CN {x,t)I{x,t). 



(17) 



Given the intensity profile through an absorbing medium, 
the intensity absorbed is given by. 



dl {x, t) = —jil (x, t) N (x, t) dx, 
which can be rewritten as 

/ (x, t) N (x, t) = (x, t) . 

fi dx 



(18) 



(19) 



Substituting Equation [TH] into equation [T71 and defining 



a 



yields. 



Hs {x, t) = -Oi-^I (a;, t) . 



(20) 



Using Equation [501 Equation [T3] can then be written 



as 



dl{x,t) _ 
dx 

d (P 

^Tlce {x,t) - -F^lce^Tlce {.x,t) . 



(21) 



which can be used to determine fully the thermal spatial 
strain profile as a function of time of the LCE. 

In previous experiments, shown in Figure[Tl the LCE is 
sandwiched between two glass slides. As such, heat will 
flow through the glass, then the surrounding air. Given 
the large volume of glass, heat will flow from the LCE 
to the glass for the duration of the experiment while the 
laser is on, and heat convection will occur across the sur- 
face of the glass-air interface. This effect will be included 
in the LCE heat equation to give a more accurate descrip- 
tion of the LCE's transient temperature profile. 



We assume that the one-dimensional heat equation ap- 
proximates the sample and the glass slide with contact 
conduction boundary conditions at the illuminated LCE- 
glass interface and convection at the other end of the 
glass slide. Therefore, Equations [T3l through [TBI will de- 
scribe the temperature for < x < Lqi and the heat 
equation for the glass substrates with the step function 
intensity described earlier will be modeled using the fol- 
lowing equations for — Lg < x < 0. 



d 

^Tglass {x,t) — ii'glass Tglass (x,t) , (22) 



d 

dux 

d 
dux 

C 



X — — L r, 



[Tglass (0,t)-TLCE (0,t)], 



-T^uss{-Lg,t), (23) 



^glass 
Tglass (X,0) =0 



(24) 
(25) 



where fcgiass is the thermal conductivity, ifgiass is the 
thermal diffusivity, and Lg is the thickness of the glass, 
while h is the convection coefficient between glass and 
air at room temperature. Note that the glass has only 
limited interactions with any substance other than the 
LCE and air. Equations through will also be used 
to find the temperature throughout the far glass slide at 
the dark interface located at Lq < x < Lq + Lg, where Lq 
is the coordinate of the LCE-glass interface and Lq + Lg 
is the glass slide's surface that is in contact with air and 
cools via convection. 

Once the temperature is known throughout the LCE, 
the strain induced by the change in temperature can be 
found. In the case of a nematic LCE, the scalar uniaxial 
order parameter is also temperature dependent. Thus, 
in the 1/e absorption length, photo-isomerization could 
dominate the photomechanical effect. However, the heat 
generated in this region can flow to the dark region, lead- 
ing to a thermal photomechanical response. 

The thermal strain, at (x,t), is defined as 



ALt 

Lo ' 



(26) 



where ALt is the change in length caused by the thermal 
processes. 

Similarly, following the previous approximations 
for photo-isomerization-induced strain, the thermally- 
induced strain can then be written as at{x,t) cx AQ 
where AQ oc AT for small changes in temperature. 
Therefore, 



atix,t) = -q[T{x,t)-T{x,0)] 



(27) 



where g is a thermal strain proportionality constant (in 
the same manner as b in Equation [5]) . 
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C. Length change 

The total strain in the decoupled approximation is just 
the sum of the two strains caused by photo-reorientation 
and thermal expansion, 



(T (x, t) = at {x, t) + CTp (x, t) . 



(28) 



The strain in Equation [21] is then used to evaluate the 
total change in length of the LCE. The change in length 
is related to the average strain, Uavg, by 



AL 



for 



AL < Ln 



(29) 



where the total length change is AL = ALp + ALt. If 
the strain were uniform, then the total change in length 
would be the strain multiplied by the initial length. In 
contrast, in our experiments the strain depends on the 
depth, so that the change in length in each segment of 
the LCE depends on its location. Therefore, if the initial 
length of each segment is Ax, then the change in length 
of the i*^ segment is ALi, which is given by. 



ALi = a{xi)Ax. 



(30) 



The total length change can be calculated by adding up 
the length of each section. 



AL = ^a{xi)Ax. 



(31) 



In the continuum limit, if the strain is sufficiently small, 
the total length can be approximated by 



AL 



u{x, t)dx. 



(32) 



Over a long time period, the LCE reaches an equilib- 
rium strain profile, accompanied by an equilibrium pop- 
ulation profile, temperature profile, and intensity profile. 
A measurement of AL as a function of time during the 
transition to equilibrium can be used to determine the 
phenomenological parameters from Equation[32]with the 
solutions obtained for Equation [55] which is used as the 
input of the numerical approximations to Equation 1211 



D. Reversibility of LCE length change 

The laser is turned off after the system reaches equilib- 
rium. Once again, the heat equation is solved but with 
no heat source; and, the initial strain profile for the re- 
laxation process will be simply the strain the instant the 
laser is turned off. 

The photo-isomerization-induced strain is treated in 
the same manner, where the population fraction of trans 
isomers along the direction of the polarization of the elec- 
tric field becomes 



dt 



N{x,t)=P{N,q-Nix,t)), 



(33) 



where the initial condition is just the population fraction 
of absorbers throughout the LCE at the time the laser is 
turned off. 

The heat equation for the glass substrates are identi- 
cal to the previous ones except that the initial condition 
is now the temperature profile at the time the laser is 
turned off. Similarly, the temperature profile of the LCE 
at laser shut-off must be used as the initial conditions 
when solving for the LCE's relaxation. The equations 
governing the cooling of the LCE are 



^TlcE (x, t) = /^LCE^TLcE (x, t) , 



dn^ 
C 



khCE 

d 
drix 
C 



■Tlce {x,t) 

[TLCE(0,t)-Tgiass (0,t)], 



Tlce {x, t) ^ 

x=Lo 



khCE 

Tlce [x, 0) = G (x) , 



(34) 



(35) 



(36) 
(37) 



where G {x) is the initial temperature profile of the LCE 
the instant the laser is turned off. 

As before, these equations can be used to calculate 
the individual strains due to the decoupled processes and 
summed. Note that during the relaxation process, we 
use the same parameters as those determined from the 
transient response after the laser was turned on. 



IV. COMPUTATIONAL RESULTS 

For the case of photothermal heating, we must take 
into account the heat generated by the laser at the input 
end as well as thermal diffusion through the sample and 
glass substrates. The photo-isomerization mechanism in 
the illuminated region is ignored due to the extremely 
short path length. Furthermore, because of the small 
strain observed in previous POD experiments, we will use 
the engineering strain approximation and will neglect the 
effects due to the constraint imposed at the LCE-glass 
interfaces. 

We solve the heat equation for the glass and LCE si- 
multaneously using the constraints due to the boundary 
and the material's constitutive relations to get a set of 
equations describing thermal-induced length changes us- 
ing the finite difference approximation. The increments 
are small enough to obtain a good approximation to our 
experimental parameters. 

The explicit finite difference method was chosen due to 
the multiple boundaries that are taken into account. The 
photo-isomerization equation and heat equation were 
solved in relatively large spatial increments, 66 /im, and 
small temporal increments, 1 ms, to account for stability 
governed by the Fourier number, which must be taken 
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into consideration for the stability of the exphcit finite 
difference approximation. Because the path length of 
the pump beam is aprroximately 5 fim for a high con- 
centration DOS-doped LCE, all of the light is absorbed 
in the first LCE increment. Thus, the change in the 
population fraction of trans isomers is negligible in the 
heat source term throughout the rest of the sample. This 
picture is consistent with experimental observations of a 
photomechanical response that is dominated by thermal 
diffusion. Using these approximations, the heat equa- 
tion was numerically solved to determine the temperature 
profile from which the length change was determined. 
The numerical models for length contraction and relax- 
ation were used as fitting functions, where q is the ther- 
mal strain proportionality constant, and is introduced as 
a phenomenological fitting parameter. 

The thermal strain, at{x,t), peaks near the surface 
of the LCE at short times, and propagates inward as 
the heat diffuses beyond the illuminated region. Photo- 
isomerization orientational hole burning leads to a de- 
crease in optical absorption, and thus allows the light to 
penetrate further into the sample. The light continues 
to burn an isomerization hole until an equilibrium inten- 
sity profile is reached, defined by the condition that the 
rate of isomerization induced by the laser matches the 
relaxation rate at each point in the sample. Since exper- 
imental evidence suggests that the photo-isomerization's 
contribution to the strain from the illuminated region is 
negligible, the step size in the numerical calculations was 
chosen to be much larger than the 1/e absorption length 
to avoid the necessity for modeling the effects of orienta- 
tional hole burning. 

Since photo-isomerization is found to occur in a thin 
surface layer of a high concentration DOS-doped LCE 
(typically less than 50 /im) , the contribution to the length 
change of the LCE from this mechanism is negligible. 
The dynamics of the photo-isomerization process, how- 
ever, determines the temporal dependence of energy ab- 
sorption; but, once this fast process reaches equilibrium, 
its net effect can be approximated as a constant heat 
source near the sample's surface. Thus, the mechanisms 
of photo-isomerization and heating are effectively decou- 
pled. 

Previous work supports the view that length contrac- 
tion due to photo-isomerization is negligible. 16| This ap- 
proximation will remain valid as long as the light is ab- 
sorbed over a thin-enough length that its contribution to 
the length change is negligible. In this regime, we can 
set N Neq, and 

a{x,t) ^ at{x,t). (38) 

The length change of the LCE can be determined from 
the angle between the plates and the length change at 
the position of the probe laser using a simple geometric 
analysis based on the geometry of similar triangles. An 
estimate of the initial phase of the interferometer at the 
start of the experiment is the largest source of experi- 
mental uncertainty. The positions of the peaks and the 



initial phase of the interferometer provides sufficient in- 
formation to determine the evolution of the length of the 
LCE during and after laser illumination. 

To determine the length change, we use the transmit- 
tance equation for a Fabry-Perot interferometer. 



/max f + Fsin2(^(AL(i) + Lo))' 

where _F is a constant, commonly known as the finesse, 
that depends on the reflectivity of the mirrors. For a 
488 nm pump laser, the separation between peaks corre- 
sponds to a displacement of 244 nm. 

The theory is used to fit the data for a range of laser 
powers of high concentration dye-doped LCEs. A few 
simplifying assumptions were made. First, the thermal 
diffusivity is assumed to be constant, independent of tem- 
perature and degree of photo-isomerization of the dyes in 
the LCE. Also assumed is that the contact conductance 
coefficient through the LCE-glass interface remains con- 
stant and is equal on both the illuminated interface and 
the opposite dark interface. The convection coefficient, 
h, is also assumed to be a constant. Finally, the beer- 
lambert coefficient is also assumed to be temperature in- 
dependent. 

The results of a numerical simulation of the tempera- 
ture over a 40 second interval is shown in Figure [SJ which 
shows the temperature as a function of time and depth in 
the LCE and glass slides. Heat flows out through all the 
surfaces of the 1 mm-thick glass substrates through con- 
vection. The substrate's large volume makes numerical 
modeling difficult due to memory and computing power 
limitations. Our approach is to apply a one-dimensional 
approximation; but, increasing the thickness out to a dis- 




FIG. 5. (Color online) The calculated temperature change 
due to heating by absorption of a 36 mW laser over the surface 
of a 400 fim LCE initially at room temperature. The middle 
peak corresponds to the position of theLCE while the broader 
wings correspond to the glass substrate. 
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tance that is larger than the actual substrate thickness 
but smaller than its two largest dimensions to take into 
account the larger volume of the substrate. While the 
thermal difFusivity of the glass is set to the literature 
value, the heat transfer coefficient is adjusted so that the 
average heat transfer through the sample/glass interface 
approximates the actual value in the presence of the large 
glass substrates. This is a good compromise because it 
takes into account heat conductance at the sample sur- 
face where the temperature is greatest and approximates 
the effects of the parts of the substrate that are furthest 
from the interface, which will not effect the sample's tem- 
perature profile significantly. In effect, this approach av- 
erages over the effects of the remote parts of the bulk 
substrate. 

Figure [5] shows the average temperature from numer- 
ical calculations and temperature probe measurements 
plotted against time. The theory and experiment agree 
well. 

The total strain approximated as the thermal strain in 
Equation [38] is used to model the total length change of 
the sample, and the numerical results are fit to the ex- 
perimental data. The length change data from measure- 
ments and the theoretical curves for several intensities 
are shown in the top of Figure [71 Note that the high op- 
tical density of the sample leads to high energy density of 
absorption, so that one must be careful to use intensities 
that are low enough to prevent damage to the LCEs. 

The values of the parameters from the literature or 
determined from fitting are shown in Table [T) Note 
that the constant of proportionality, q, for a dye-doped 
LCE is about an order of magnitude larger than it is for 
DRl-doped PMMA (approximately2.5 x 10-* K"^) as 
reported by Xiang and co-workers. [2g This difference is 
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FIG. 6. (Color online) The change in temperature for an ex- 
perimental run lasting 40 seconds (points) and the numerical 
calculation (curves) using the parameters shown in Table [T] 
The time scale is logarithmic. 
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FIG. 7. (Color online) (a) The length contraction of an LCE 
for a series of incident laser powers as a function of time after 
the laser is turned on, and (b) the length relaxation when the 
laser is turned off after being illuminated for 5 minutes. The 
theory curves use the parameters shown in Table [I] 



primarily due to the leveraged effect of a photothermally- 
induced decrease in orientational order of the mesogens. 

The parameters from the fits of the LCE length as a 
function of time in the presence of light were either de- 
termined separately by experimentation, found in litera- 
ture, or determined by fits using the numerical calcula- 
tion. The data and the fits for the illuminated response 
are shown in the top graph of Figure [7| The bottom 
graph in Figure [7| shows the theoretical curves and data 
of length relaxation to the initial state as a function of 
time after the light beam is turned off. The theoreti- 
cal curves showing relaxation use the same parameters 
as were determined from the light-induced studies. Also 
note that the set of parameters are the same for all curves 
at the various light intensities. 

The calculation of length relaxation over time after 
the laser is turned off spans 40 seconds of time. Since 
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TABLE I. LCE Parameters 



Constant Value Units Source 
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our model does not take into account the changes in 
the long-time boundary condition, the theoretical predic- 
tions for the total length contraction arc always slightly 
smaller than the experimental values. Note that assump- 
tions about the geometry of the glass and LCE leads to 
a small deviation between the experimental data and the 
theory at longer times. However, the predicted trends 
are consistent with the data. 

To elaborate, with one set of parameters, our theoret- 
ical model of photo-isomerization decoupled from pho- 
tothcrmal heating accurately predicts both the length in- 
crease as a function of time of an LCE for a wide range of 
intensities, as well as the relaxation process (Note that q 
is a function of the initial nematic order and can change 
from sample to sample). Thus, the parameters are in- 
trinsic to the system, and their dependence on the laser 
power is observed to be negligible. This self consistency 
of a broad set of data that is described by one set of pa- 
rameters that agree with literature values suggest that 
our model of the mechanisms of the photomechanical re- 
sponse is valid. 
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